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Abstract: 

In this study we explore the LHC’s Run II potential to the discovery of heavy Majorana 
neutrinos, with luminosities between 30 and 3000 fb _1 in the l±l±j j final state. Given that 
there exist many models for neutrino mass generation, even within the Type I seesaw frame¬ 
work, we use a simplified model approach and study two simple extensions to the Standard 
Model, one with a single heavy Majorana neutrino, singlet under the Standard Model gauge 
group, and a limiting case of the left-right symmetric model. We then extend the analysis 
to a future hadron collider running at 100 TeV center of mass energies. This extrapolation 
in energy allows us to study the relative importance of the resonant production versus gauge 
boson fusion processes in the study of Majorana neutrinos at hadron colliders. We analyze 
and propose different search strategies designed to maximize the discovery potential in either 
the resonant production or the gauge boson fusion modes. 
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1 Introduction 

There is convincing evidence for the existence of three active neutrino species [1, 2], at least 
two are massive. Their mass pattern has been narrowed down to a normal or an inverted hi¬ 
erarchy or degenerate [3-8]; however, the absolute mass scale remains unknown. In addition, 
we are now capable of measuring CP violation in the lepton sector, given that the last mixing 
angle, $ 13 , has been measured by several reactor experiments [9-11] and the T2K accelerator 
experiment [12]. The experimental status of neutrino physics suggests that the evidence for 
neutrino masses represents a clear motivation for new physics beyond the Standard Model 
(SM). In principle this is due to the fact that within the SM, neutrinos are massless and one 
may incorporate new degrees of freedom or effective interactions to generate Dirac or Ma¬ 
jorana masses. The latter is an interesting possibility since a Majorana mass term violates 
lepton number by two units. While a neutrinoless double beta decay (Oz^/3/3) signal would be 
a breakthrough giving us knowledge on the nature of neutrinos and with additional assump¬ 
tions the scale of the active neutrino mass matrix, the observation of lepton flavor violating 
processes at the LHC can shake the foundations of the SM with additional mechanisms and 
with luck we may directly find the new particles associated with these interactions. 
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In the minimal Type I seesaw mechanism [13], the SM is extended with a single Majorana 
fermion, singlet under the SM gauge group. Within this framework, the Majorana fermion 
carries lepton number and couples to left-handed leptons through the Higgs field. In addition, 
a Majorana mass term for the Majorana fermion can be implemented consistent with the 
gauge symmetries of the SM. The Majorana nature of this fermion and its mixing with the SM 
active neutrinos may lead to interesting testable interactions at the LHC that violate lepton 
number such as the production of same-sign leptons in association with jets. This final state 
can be achieved through resonant production of a Majorana neutrino and through gauge boson 
fusion. There has been a number of theoretical studies aimed at examining the sensitivity 
that the LHC has to heavy Majorana neutrino masses and to set limits on the couplings to 
leptons that arise from mixing with the active SM neutrinos [14-21], In particular, a 14 TeV 
LHC with 100 fb -1 may have sensitivity to resonant production of Majorana neutrinos with 
masses up to ~ 400 GeV using the and n ± e ± channels. However, the VFy fusion can 

reach masses up to 530 GeV in the channel [22], Furthermore, in the regime where the 

Majorana Neutrino has mass below the mass of the SM W boson, the authors in [23] have 
shown that the LHC can be sensitive to mixing angles in the range 10~ 4 —1CP 3 and Majorana 
neutrino masses below 15 GeV. This is done by looking for a prompt lepton in association 
with a displaced lepton jet. For heavy Majorana masses, no significant reach can be achieved 
with the gauge bosons fusion process. These channels are collider analogue to (Oz//3/3) which 
are not available there. Moreover, the ee channel may be sensitive to electron-like Majorana 
neutrinos with masses between 150 — 1400 GeV and higher LHC energies [24]. 

Alternatively, one can extend the SM to be left-right symmetric; partnering right-handed 
charged leptons with new right-handed neutrinos into an SU(2 )r doublet. Within this class 
of models, the gauge couplings are left-right symmetric and a new charged gauge boson, Wr , 
connects right-handed charge currents [25-29]. The phenomenology of left-right symmetric 
models is very rich with the appearance of new scalar degrees of freedom and lepton violating 
process such as the production of same-sign leptons. The latter can be induced from the 
production of a right-handed charged gauge boson. A number of theoretical studies have 
examined the sensitivity that the LHC has to elements of the right-handed lepton mixing 
matrix, the mass of Wr and Majorana neutrino masses [30-41]. 

With the LHC running, both CMS and ATLAS have searched for same-sign dilepton 
final states using light- or 6-jets [42-44]. In addition, a recent analysis by CMS [45] is used 
to search for heavy Majorana neutrinos in the same-sign dimuon channel with an integrated 
luminosity of 19.7 fb _1 and 8 TeV center of mass energies with a sensitivity to a Majorana 
mass of 500 GeV. Dedicated searches to probe same-sign leptons within the framework of 
left-right symmetric models have also been carried out. In particular, both CMS and ATLAS 
have set limits on heavy Majorana neutrino production assuming M\y R > M]\r R and identical 
quark and neutrino mixing matrices for the left- and right-handed interactions [46, 47]. The 
exclusion regions extending to M\y R ~ 3 TeV in the (M\y R , Mr^ It ) plane. 

Our work compliments and extends the various analyses discussed above. In the present 
analysis we explore the LHC’s Run II potential to discover heavy Majorana neutrinos with 
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luminosities between 30 and 3000 fb _1 the l±l±j j final state. Given that there exist many 
models for neutrino mass generation, even within the Type I seesaw framework, we use a 
simplified model approach and study two simple extensions to the Standard Model, one with 
a single heavy Majorana neutrino, singlet under the Standard Model gauge group, and a 
limiting case of the left-right symmetric model. We would like to emphasize that simplified 
models are not complete models. They are constructed to highlight specific points which 
in our case means LHC signals, and they can be mapped into more realistic models. They 
can also be viewed as truncations of more complete models. For the Type I Seesaw, a 
more natural framework would include at least two heavy right-handed neutrinos to explain 
the mass differences between the active neutrinos. However, in the mass eigenstate basis 
of these new heavy Majorana neutrinos the mixing angles and the masses can be treated 
as independent variables. Furthermore, more details regarding the model are necessary to 
determine the correlation between the new parameters and the mass differences of the active 
neutrinos. We will be conservative and assume that only the lightest of the new heavy sates 
are within reach of the LHC or future circular collider. Incorporating a second and third 
heavy Majorana neutrino simply leads to a duplication of our results; thus we ignore the 
heavier degrees of freedom and also the mass relation between active neutrinos. Thus, we 
can focus on the lightest heavy Majorana neutrino without loss of generality. Furthermore, 
the usual Type I seesaw model has very small mixing of light and heavy neutrinos. This can 
be circumvented by arranging structures in the heavy neutrino mass matrix. The details are 
highly model dependent; albeit important for fitting low energy neutrino data, see e.g. [48, 49]. 
In the simplified model this is captured by allowing the couplings of the heavy neutrino to the 
charged leptons to be free parameters. We then extend the analysis to a future hadron collider 
running at 100 TeV center of mass energies. This extrapolation in energy allows us to study 
the relative importance of the resonant production versus gauge boson fusion processes in the 
study of Majorana neutrinos at hadron colliders. We analyze and propose different search 
strategies designed to maximize the discovery potential in either the resonant production or 
the gauge boson fusion modes. Our work is strictly phenomenological. We check consistency 
with the strongest model independent constraints that arise from low energy interactions such 
as rare decays of the muon, unitarity of the PMNS matrix, and the null evidence for (Ov/3/3). 
A well known caveat to be added here is that other than unitarity tests, low energy constraints 
are obtained with the assumption that other new physics do not contribute significantly. 

The summary of our study is as follows: In Section 2 we introduce the models and 
interpret current low energy and collider constraints. In Section 3 we calculate the leading 
SM backgrounds and propose various search strategies designed to maximize the discovery 
potential in either the resonant production or the gauge boson fusion modes of a heavy 
Majorana neutrino. In Section 4 we provide concluding remarks. 
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2 Two Simplified Models 


To study same-sign leptons production in association with jets at the LHC run II and future 
higher energy colliders, we examine two simplified models that give rise to this final state 
topology. The first one consists only of SM gauge interactions and the only added new 
degree of freedom is the singlet heavy Majorana Nr fermion which mixes with the 3 active 
SM neutrinos. The mass of Nr and its mixing with the active neutrinos are taken to be free 
parameters. The second model implements an additional SU(2) gauge symmetry under which 
the right-handed leptons and quarks are charged unlike the left-handed chirality. Nr now 
partners with the right-handed charged leptons. This can be part of a left-right symmetric 
model that can give rise to realistic masses for the active neutrinos. The details are not 
essential for our purpose. In what follows, we review different aspect of these two theoretical 
approaches and indicate how they can be related to more complete models. 

2.1 Single Majorana Fermion Extension to the Standard Model 

The following framework consists of an extension to the SM with a single Majorana fermion, 
singlet under the SM gauge interactions. Within this analysis, the Majorana fermion, Nr, 
mixes with the three active neutrinos. We are aware that to get realistic active neutrino 
masses, two or more Majorana neutrinos have to be added. In the spirit of the simplified 
model we are considering we can assume that these additional states are either much heavier 
or have negligible mixings. Following the analysis detailed in [17] where in terms of mass 
eigenstates, the gauge interaction Lagrangian is given by: 

c = -~k w * - jsK IE 

* \i=em=l J * \l=e ) 

< 2 '« 

where U and V provide a relationship between the flavor basis, v' lL and the mass basis 
parametrized by v m L such that 


3 

u lL = ^2 Ui m v m + ViRfN^, (2-2) 

m— 1 

with UU^ + VV^ = 1. We refer the reader to [19] for more details regarding the mixing 
formalism. In our simplified approach, we parametrize the charged current interactions by 
introducing the parameters e/, where 

e; = Vin- (2.3) 

For Type I seesaw models the mass of the Majorana neutrino is commonly assumed to be 
Mn r > 10 14 GeV to obtain active neutrino masses below 2 eV [50]. However, there exist 
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well motivated models that can lower the seesaw scale and/or have a sizable coupling to the 
SM particles. Notable examples are structures or cancelations in the heavy neutrino sector 
that can lead to Mn r ~ O (TeV) and larger mixings [49]. Approximately conserved lepton 
number has also been invoked [51, 52], The inverse seesaw mechanism can also yield TeV 
scale Nr masses with mixings at a few percent [53]. In our analysis we will take a model 
independent approach and vary the e couplings freely to investigate the range that the LHC 
can probe. 

(b) 


u £+ £'+ j l 



u £+ £'+ j l 



Figure 1. Feynman diagrams for the ^-channel contribution to same-sign lepton production. 


The signal that will be probed in this study is given by 

pp^ Z± Z* + j j, (2.4) 

where l, l 1 denote either an electron or muon. One important aspect of this reaction is that 
it contains no missing transverse energy at the parton level. Most of the missing energy will 
appear after a full detector simulation has been carried out, mostly due to the misidentification 
of jets. At the parton level there are two channels to consider: 

1. s-channel ud annihilation depicted in Figure 1: 

u +d —y W* — >l + + N r —> l + l ,+ + W, (2.5) 

with the W decaying into light jets. The Majorana neutrino, Nr can be on- or off-shell 
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depending on the mass, Mn r . In this channel, the W from the decay of Nr can allow 
us to control the large SM backgrounds since one may reconstruct the W mass using 
the two light jets for not too large Mn r . 

2. t-channel or WW fusion process depicted in Figure 2: 

u + u —» d + l + l' + + d. (2.6) 

In this process, the Majorana neutrino is always off-shell. If the initial quarks have the 
same color, a contribution to the amplitude will arise by interchanging the two forward 
outgoing jets. 


“(pi) 




(d) M d 


Figure 2. Feynman diagrams for the T-channel contribution to same-sign lepton production. Dia¬ 
grams (c) and (d) correspond to contributions where the final state patrons have the same color. 

2.2 Gauge bosons beyond the Standard Model 

The signal topology described above is not unique to models that generate mixing between a 
heavy ”right-handed” Majorana neutrino and the three active neutrino species. In fact, mod¬ 
els that go beyond the SM and that include additional gauge degrees of freedom can give rise to 
same-sign lepton final states. Left-right symmetric models are primarily motivated to address 
the origin of parity violation in low energy physics. It assumes that the Lagrangian is left-right 
symmetric, with the observed asymmetry, such as in /3-decay, arising from a vacuum expec¬ 
tation not invariant under parity transformations. The SU(2)lxSU(2)rxU(1)b-lxPl.o-R) at 
low energies, breaks down to the SM gauge group with new particles mediating new inter¬ 
actions arising at higher energies [25-29], in particular, right-handed charged currents, Wj^. 
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Below we summarize essential ingredients pertaining to left-right symmetric models in light 
of a same-sign lepton final state at the LHC and refer the reader to [54] for a full review of 
the theory. 

In our analysis we assume that new charged currents are associated with a gauged SU(2 )r 
symmetry and introduce three heavy Majorana fermions, N ^' T , that together with the right- 
handed charged leptons transform under the fundamental representation of the new gauge 
symmetry. The matter content can be written as 

= t)« : ( 1,2 ’ 1/3 )’ 


Ql = 


: (2,1,1/3), 


Qr 


Ll 



■ ( 2 , 1 ,- 1 ), 


Lr 



( 1 ) 2 , — 1 ). 


(2.7) 


Here by symmetry we have three right-handed neutrinos N l . Within this framework, the 
minimal setup consists of one scalar bidoublet and two complex scalar triplets that are given 
by 


4 > = 



Al,r 


61,r/V2 \ 

$L,R S + L,r/V2) 


( 2 . 8 ) 


This is required to break the gauge symmetries down to U(l) em . The scalar phenomenology 
is very rich and recent studies have been devoted to the possibility that a doubly charged 
scalar can be detected at the LHC and future more energetic hadron colliders [55-57]. In 
addition, production and decays of charged scalars could help determine the leptonic right- 
handed mixing matrix at the LHC [40]. In principle, the heavy gauge bosons will mix with 
the SM gauge sector, but this mixing is naturally small due the V-A nature of the charged 
current observed at low energies. In our analysis we explore the limit where the mixing is 
negligible. 

In the minimal setup, lepton masses are due to the following Yukawa interactions 


Cy 


= L l Lr + i (LlCia 2 Y AL A L L L ) 

+ 2 ( LrC ia 2 Y A R ArLr) . 


(2.9) 


In the absence of spontaneous CP violation we obtain a Dirac mass term mixing left- and 
right-handed leptons given by 

Mu = vi Y<j> + Y$v 2 , (2-10) 

where v\ and v 2 are the vacuum expectation values (vev) of the Higgs bidoublet, and cj) 
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In addition both left- and right-handed neutrinos acquire a Majorana mass given by 


Mu 

m n r 


Y Al v l - Mp——M D , 

-MNr 

h a r vr, 


( 2 . 11 ) 


with vl = (<5l) and vr = {Sr). In contrast to the Type I seesaw models vl need not be zero 
at the tree level, although it is expected to be small compared to the Fermi scale. The details 
are model dependent. The diagonalization of the mass matrices proceeds as in the Type I 
seesaw, leading to charged current interactions that in the mass eigenstate basis are given by 

Ccc = (vlVlWlIl + N 1 r VrWrIr) + H.c, (2.12) 

where Vl r are the 3x3 left- and right-handed mixing matrices respectively. In the above 
Lagrangian, Vl and Vr have been obtained in the limit where Md/Mn r <C 1. In this 
limit Vl = E/pmns an d parametrizes the mixing of light neutrinos through charged current 
interactions leading to neutrino oscillations, while Vr parametrizes new phenomena that can 
be probed at the LHC and through various lepton flavor violating processes [40, 41]. 

The goal of this study is to probe Vr, in particular the signal topologies in Figures 1 
and 2 with W —>• Wr at the LHC with 13 TeV center of mass energies and at a future 100 
TeV Collider. To accomplish this, we focus only on the structure given in Equation (2.12) to 
parametrize the charged currents in the lepton sector and assume that Wr couples at leading 
order to quarks with a hierarchical structure similar to that of the SM CKM matrix. This is 
done to suppress the bounds on Mw R that arise from measurements of CP violating effects 
of the Kl — Kr mass difference and BB s meson parameters [58-60]. However we direct the 
reader to [61-63] where a general study on right-handed quark mixings can be found. 

2.3 Constraints 

In the singlet Majorana extension to the SM introduced Section 2.1, the Majorana neutrino 
mixes with active neutrinos with a strength proportional to e/ for l = e,fi, r. As a result, this 
framework is sensitive to model independent constraints that arise from lepton unitarity [64, 
65]. These are 

|e e | 2 < 0.012, |e M | 2 < 0.0096, |e T | 2 < 0.016. (2.13) 

In addition, indirect constraints that are model dependent arise from rare decays of the muon 
and muon properties, such as /r —>• and the muon anomalous magnetic moment, a^. The 
validity of these constraints assumes that Nr is the only source of new physics entering the 
calculation of the matrix elements, such that there is no cancellations arising from unknown 
physics. The Feynman diagram is depicted in Figure 3 and we use the unitary gauge to avoid 
the need to implement the coupling of Nr to the would-be goldstone modes. The amplitude 
is given by 


T = A^{q)u e {p- g)[«V,^(l + 7 5 )]?v(p), 


(2.14) 
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Figure 3. Contribution to the g —> ey decay rate in the presence of a heavy Majorana neutrino. In 
left-right symmetric models, the the SM W gauge boson can be replaced by Wr. 


where £ /J ' denotes the photon polarization vector with momentum q. Charged lepton self¬ 
energy diagrams and the active neutrino diagrams that are necessary to cancel divergences 
are not shown. The Majorana neutrino contribution is given by 


A = 


i g 2 e rrin 
32vr 2 8 M 2 y 


Cfi^eF 




(2.15) 


where F(a) can be determined from 


F(a) 


13 

¥ 

13 

¥ 


a(l — 6a + 3 a 2 + 2a 3 — 6a 2 In a) 

/(«) 


(2.16) 


where a = /M'w ■ The contributions from the active neutrinos can be obtained by taking 

the limit Mn r —>• 0. These contributions can be subtracted out by using the unitarity of the 
neutrino mixing matrix. With the above two equations, the branching ratio for g —>• ey is 
given by 


BR(n —>• ej) 


3 a 

3 ¥^ 

3 a 

s |e 


i/(»)l : 


(2.17) 


where the second equality is valid in the limit of very heavy Nr compare to the gauge boson. 

For the case of right-handed gauge interactions we use the unitarity of the right-handed 
neutrinos mixing matrix and take the limit where two of the three N l R decouple from the low 
energy theory to obtain 


BR{[i —»• ey) 


3a 

3 ¥ 


Mw_\ 

m w r ) 


4 


\V R RiV* Rel 


2 |2-/(a)| 


2 


(2.18) 
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where M\y denotes the mass of the SM charged gauge boson and {M\y /M\y n ) A arises due to 
normalization. The current experimental bound is BR(/i —>• ey) < 5 x 10~ 13 [66]. 

Similarly, the contribution from Nr to the anomalous magnetic moment of the muon, a p 
if given by 


A a p 


5 15 — 74a + 114a 2 — 54a 3 — a 4 — 12a 2 (2 — 3a) In a 

12 + 12(1 - a) 4 

^c, a —y 0 
|c, a —>• oo 


(2.19) 


where c = |e^| 2 . The first limit in the above equation is consistent with the results 

appearing in [67]. The experimental limit for the deviation is a)) xp — a tAl = 287(63) (49) x 10 -11 . 
A few comments are in order: 


1. Neutrinos, heavy or otherwise, contribute negatively to a p and will not explain the 
alleged discrepancy with the SM. 

2. The calculation of A a p does not make use of the unitarity condition of the lepton mixing 
matrix. 


Neutrinoless double beta decay can be a sensitive probe for the absolute neutrino mass 
scale and/or the mixings of heavy neutrinos with the active ones. The caveat here is that the 
theoretical uncertainties involved in the calculation of the nuclear matrix elements are difficult 
to quantify. To compound this when heavy neutrino intermediate states are involved one is 
dealing with very short distances between the two neutrons involved, and the uncertainties 
are even larger. With this in mind, the best limit on this reaction comes from the decay of 
76 Ge which gives an upper bound on the Major ana active neutrino mass of rn v < 300 — 600 
meV [68]. This bound translates to a bound on the light-heavy mixing element V e jy = e e > 
Section 2.1, and it is given by 

< 5 x 10“ 8 GeV" 1 . (2.20) 

M N R 

Within left-right symmetric models, Section 2.2, a similar bound can be obtained on the 
mixing matrix elements between the right-handed leptons, V e i, where i runs over the three 
heavy right-handed neutrinos, A/r e,7l ’ r . In particular, one can constrain 


e N R 


3 


T, v « 


m p Myy_ 

M Nr M W r 


( 2 . 21 ) 


[41] using the current experimental limit of |ctv | < 2 x 10 8 [69]. 
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In addition to low energy direct and indirect constraints, collider searches can be used to 
probe the single Majorana extension of the SM, in particular flavor violating final states. The 
most recent analysis by the CMS collaboration targets the same-sign dimuon channel with 
yfs = 8 TeV center of mass energies and 19.7 fb _1 of integrated luminosity [45]. The limits 
obtained are on the V^n mixing element and result on the following upper bounds: | | 2 < 

0.00470, 0.0123,0.583 for Majorana neutrinos masses of 90, 200 and 500 GeV respectively. 
However, for a 200 and 300 GeV Majorana mass the unitarity bound given in Eq. (2.13) is 
stronger, making the latter the primary constraint limiting the value of that we use in 
our analysis. This analysis equally constrains the left-right symmetric model discussed in the 
previous section. In addition, both CMS and ATLAS have reported limits on heavy Majorana 
neutrino production in the context of the left-right symmetric models [46, 47]. The limits 
are for heavy Majorana production with M\y n > Mn r and identical left- and right-handed 
mixing matrices. The region excluded extends to M\y n ~ 3 TeV, with a region defined in the 
M\v R — Mn r plane. 


3 Same-sign leptons at the LHC and beyond 
3.1 Background 

Within the SM, lepton number violating processes are absent at zero temperature. However, 
there exist certain SM processes that can give rise to a same-sign lepton final state in as¬ 
sociation with jets. The two leading backgrounds discussed below are simulated at leading 
order using MadGraph 5 [70]. We implement PYTHIA [71] for the parton showering and 
hadronization and a fast detector simulation is carried out using Delphes 3 [72]. The de¬ 
tector simulator is used for jet clustering, track reconstruction and lepton identification and 
isolation. Light jets are reconstructed using the anti-fey algorithm with parameter A R = 0.5. 
Furthermore, we calculate the lepton isolation criteria, //, by summing the px of all tracks 
within a cone of size A R = 0.5 around the leptons and require that 


T Ei^lPT, 

h =- 

PT,l 


< 0 . 10 . 


(3.1) 


The two leading background processes at the LHC are the following: 

1. p p W + W + j j with low missing transverse energy, 


This SM process is an irreducible background and dominates when the amount of miss¬ 
ing transverse energy is small and comparable to the signal. The parameter space for 
heavy Majorana neutrinos within the SM was studied in [17] using 14 TeV center of 
mass energies and a luminosity of 100 fb . The authors identified the following series 
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of cuts that would suppress this background and enhance the significance of a pEp±j j: 

Pt,h > 10 GeV, < 2.5 

pr.j > 15 GeV, rjj < 3.0 

A Rf] in > 0.5, E t < 25 GeV. (3.2) 

However, a search performed by the CMS collaboration for same-sign muons and elec¬ 
trons using 4.9 fb -1 of data at y/s = 7 TeV [73] determined that the leading background 
for same-sign leptons in association with jets is from the QCD multijet component with 
misidentify jets as leptons. Consequently, the impact of the QCD background should 
be analyzed further. 

2. QCD multijet background. 

This background is mostly due to p p —>• 4 j with two jets misidentified as leptons. 
In addition, leptons may arise from decays of heavy flavored jets. Although there are 
studies by ATLAS on muon fake rates [74], the electron fake rate does not have a ref¬ 
erence value. However, from the studies in [73] and [74], the electron fake rate can 
be anywhere between 10 -4 — 1CP 5 . In what follows we refer to the fake rate as 
and emphasize that the fake rate highly depends on the detector and energy, making it 
difficult to determine with a fast detector simulator such as Delphes. In our analysis, we 
will assume a given fake rate and determine the background suppression as a function of 
the misidentification efficiency. We extract an upper bound on this rate by comparing 
our simulated background after applying the following kinematic cuts that we find to 
significantly enhance our signal over the irreducible background described above: 

Pt.jj, > 10 GeV, < 2.5 

pr.j > 20 GeV, rjj < 5.0 

m h,L > 10 GeV, Et < 40 GeV, (3-3) 

where l\ and h denote the leading an sub leading leptons, to the number of simulated 
events using the latest search for pEpE by the CMS collaboration [73] with 19.7 fb” 1 at 
yfs = 8 TeV. The result is a muon misidentification rate (fake rate) of ~ 2 x 10~ 5 . 
We use a generous upper bound for the electron misidentification rate of 10~ 5 . 

In Tables 1 and 2 we show the generated backgrounds at parton level and after a full fast 
detector simulation for yfs = 13,100 TeV center of mass energies respectively with 3000 fb _1 
of integrated luminosity. The detector level results are shown applying the cuts discussed in 
Equation (3.3) and the requirement that the leading two jets arise from the decay of a SM 
W gauge boson. The latter is used to enhance the sensitivity to the single heavy Majorana 
neutrino extension to the SM discussed in Section 2.1. 
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e + e + uu jj 

e + R + uv jj 

r+r+vv jj 

4 3 • e 'Ui 

Parton level 

cr (fb) 

0.3927 

0.7849 

0.3927 

1.79 


N 

1178 

2354 

1178 

5370 

Detector level + (3.3) 

a (fb) 

0.1187 

0.2674 

0.1187 

0.471 


N 

356 

802 

356 

1410 

Detector level + (3.3) 

+ (M\y — 10.) < m 7 - 1l7 - 2 < {M\y + 10.) 

cr (fb) 

0.0094 

0.0202 

0.0119 

0.0643 


N 

28.3 

60.7 

35.6 

193 

Detector level + (3.3) 

+{M W — 5.) < < {M\y + 5.) 

a (fb) 

0.0047 

0.0100 

0.0059 

0.0323 


N 

14.18 

30 

17.6 

96.9 


Table 1. SM backgrounds at 13 TeV. The event number is shown at a luminosity of 3000 fb 


-l 




e^e^vv jj 

e + R + vv jj 

H + n + vv jj 

4j • e Ul 

Parton level 

cr (fb) 

3.686 

7.370 

3.694 

44.2 


N 

11058 

22110 

11082 

133000 

Detector level + (3.3) 

a (fb) 

0.962 

2.148 

1.189 

11.6 


N 

2887 

6444 

3568 

35000 

Detector level + (3.3) 

+(M\v — 10.) < m 7l , 7 - 2 < (M\y + 10.) 

a (fb) 

0.0597 

0.1388 

0.0784 

1.59 


N 

179 

416 

235 

4770 

Detector level + (3.3) 

+(Mw — 5.) < m,-, , 7 - 2 < {Mw + 5.) 

a (fb) 

0.0291 

0.0692 

0.0395 

0.798 


N 

87.2 

207.6 

118.4 

2390 


Table 2. SM backgrounds at 100 TeV. The event number is shown at a luminosity of 3000 fb 


3.2 SM + Nr: S-channel collider reach at 13 and 100 TeV 



Figure 4. The S- and T -channel contribution using default parton-level cuts implemented in 
MadGaph 5 at 13 (left) and 100 (right) TeV using a flavor universal mixing angle of e = 0.1. 

The presence of heavy Majorana neutrinos in either a structure only consisting of SM 
gauge symmetries or an extension of the SM with an additional non-abelian group structure 
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will contribute to the process p p —>• jj in two topological classes: S- and T-channels. 

The S'-channel topology is depicted in Figure 1 and it is dominated by quark anti-quark 
fusion, in particular, the first generation. The T-channel topology is depicted in Figure 2 
and has contributions from both quark-quark and quark-antiquark fusion. In the massless 
limit, for quarks, these two contributions have different helicity amplitudes. In this section 
we examine the sensitivity that Run II at the LHC and a future 100 TeV collider will have 
to the model discussed in Sections 2.1. For the left-right symmetric model, a combined S- 
and T-channel analysis is performed, since within this theoretical framework, the mass of 
the SU(2 )r gauge boson is also a free parameter. The simulation of the signal at the parton 
level is carried out using MadGraph 5 [70]. In Figures 4(a) and 4(b) we show the S- and T- 
channel contributions to the same-sign electron cross section for the LHC running at y/s = 14 
TeV (left) and y/s = 100 TeV (right) within the framework a single heavy Majorana neutrino 
extension of the SM using a coupling of V e N = 0.1. From the plots we can observe that the 
T-channel contribution to the production cross section starts dominating for masses, M^ R , 
approximately above 1 TeV. 

In the above section we introduced the two most dominant backgrounds that must be 
suppressed to enhance the sensitivity that a hadron collider requires to probe a model with 
a lepton number violating same-sign lepton final state. Below we examine and introduce a 
search strategy aimed at extracting a statistical significant signal for Majorana masses below 
1 TeV. The constraints discussed in Section 2.3, in particular p —>• ey, narrow down the 
relevant final states to e±e±j j or p±p±j j, since in either case the cross section depends 
only on e e or e^, and we can use either coupling to suppress the p —» ey constraint. 

The simulation of the signal at the parton level is carried out using MadGraph 5 [70] with 
model files generated with FeynRules [75]. We use two input parameters in our simulation, 
the mass of the heavy Majorana neutrino, Mj\r R , and a universal mixing to the active neutrino 
species, e e = = e = 0.1. We separately generate e ± e ± j j and p±p±j j final states. In 

addition we focus on Majorana masses above 100 GeV to avoid constraints from LEP [76-78] 
that restrict || 2 to values below 10 -4 — 10~ 5 for Majorana masses between 5 — 80 GeV. The 
partial widths of Nr are thus given by 


^N-^l+W~ ~ ^N-^l~W+ ~ ^N^uZ 


„2 2 

■J_ _j/ 2 

64vr A Nr 


Mw 


+ 1 


(3.4) 


We begin the analysis by applying the cuts used to treat the backgrounds described in 
Equation (3.3). However, the full detector simulation misidentifies a number of leptons as 
jets. The survival rate for our signal is shown in Figure 5(a) for a 13 TeV collider and in 
Figure 5(b) for 100 TeV. A clear signature of this model is the W that arises from the decay 
of a Majorana neutrino. Since we are looking at the hadronic decay mode of the gauge boson, 
the W can be tagged using the invariant mass of two jets in the event, in particular the 
leading two jets. In Figure 6 we show the invariant mass for different masses of the Majorana 
neutrino, Mn r . 
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Figure 5. Survival rate for a same-sign lepton final state with at least two jets for 13 TeV (left) and 
100 TeV (right) center of mass energies. 



invariant mass of jij2 [GeV] 


Figure 6. Invariant mass distribution of the two leading jets in our eejj signal for various heavy 
Majorana neutrino masses, Mn r . 


To analyze the reach at the LHC and in a future 100 TeV collider to the S-channel signal 
of the simplified model, we define a significance variable by 


S 

VsTb 


(3.5) 


where S denotes the number of signal events and B the number of background events. Since 
we are using an inclusive final state (in charge) with same-sign leptons, the production cross- 
section is given by 


a oc e 2 [T(u, d -> W*)BR(N -> l + + jj') + T(u, d -> W*)BR(N -> l~ + jj')\ 

= e 2 [T(u, d —>• W*) + T(u, d -> W*)\ BR(N -> l* + jj'). (3.6) 

Given that Nr is Majorana, BR(Nr —» / + + jj') = BR(Nr —>• l~jj') and given the fact 
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that they are independent of the momentum transfer, we can treat e 2 BR(N —>• + jj r ) as 
a scaling factor to the total amplitude. Therefore, we choose to carry out the simulation for 
a fixed value of e = 0.1 for all three lepton flavors. Thus 


e 2 BR(N e ± + jj') e=Vei=Vlii =v Ti = \V e4 \ 2 BR(N ->■ e ± + jj'), (3.7) 


and since BR{N —>• + jj') 

write 


= 2/27 for the range of masses that we are considering, we can 
^e 2 = \V eA \ 2 BR(N ^e ± + jj'). (3.8) 


Similarly for the nfijj final states. In our analysis we include systematic and statistical 
uncertainties and model the distribution as Gaussian with an statistical error given by yN, 
for N number of events. We find the 95% confidence level fit to e using a x variable defined 
by 


X = 


e 2 S 0 


V^So + Bfr + Bsm + cr S y S 


(3.9) 


where e 


e/0.1, S 0 denotes the number of signal events with the benchmark value e = 0.1 




(a) (b) 


Figure 7. The 95% expected limit at 13 TeV and 3000 fb _1 of integrated luminosity. On the left, the 
limits are shown without the QCD contribution while on the right, varying fake rates are introduced 
for the e ± e ± jj final state. 


and Bsm,Bfr denote the SM irreducible and QCD backgrounds respectively. The latter 
is scaled with the fake rate described earlier in the section. The systematic uncertainty is 
parametrized by a sys and we assume it to be negligible compared to the statistical compo¬ 
nent. In Figures 7(a) we plot the 95% C.L region of parameter space against the irreducible 
background using a center of mass energy of 13 TeV and 3000 fb -1 of data. In Figure 7(b) 
we show the effect the QCD background with different fake rates in the eejj channel. In 
Figures 8(a) and 8(b) we show the results of the simulation for a 100 TeV collider with 3000 
fb -1 of data. It is worth mentioning that our results improve by a non-negligible amount 
if instead of doing our fit to the the whole range of Majorana masses, we carry out a fit to 
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(a) (b) 

Figure 8. The 95% expected limit at 100 TeV and 3000 fb” 1 of integrated luminosity. On the left, the 
limits are shown without the QCD contribution while on the right, varying fake rates are introduced 
for the e ± e ± jj final state. 


bins of the Majorana mass parameter, Mjy R . In particular, we define a bin size of 10 GeV 
and fit our simulated data to obtain the 95% C.L regions; these are depicted in Figures 9(a) 
and 9(b) together with the results shown in the previous two figures for comparison. It is 
worth pointing out then that our analysis shows that a 13 TeV machine will have sensitivity 
to couplings approximately an order of magnitude below those probed by the latest CMS 
same-sign muon search [45]. 




Figure 9. The 95% expected limits with a fit to 10 GeV bins of the invariant mass, Mijj. The plot 
on the left corresponds to 13 TeV center of mass energies with varying fake rates while the plot on 
the right is a 100 TeV machine. For comparison the unbinned fit is shown in red. 


3.3 SM + Nr: T-channel collider reach at 13 and 100 TeV 

We saw at the beginning of Section 3.2 that the 5-channel contribution to the production 
cross section decreases significantly for Majorana masses above ~ 1 TeV and the T-channel 
contribution begins to dominate. This effect can be easily understood given that for quark- 
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antiquark invariant masses below the mass of the heavy neutrino, the pp —>• l + Nr is closed 
while the pp —>• l + l + W~ is open; that is 

« > M N iV O’ (pp -> l + l + jj) dominant = ° {PP ^ l + Nr) BR(Nr -> /+ jj) 

S < M Nri ° (PP^ l +l+ Jj) dominant = (?>P l+l+ W ~) BR{W jj). (3.10) 

Although the phase space contribution from a 3-particle final state is larger than a 4-particle 
final state (T-channel), the S'-channel amplitude is greatly affected due to the PDF suppres¬ 
sion at high momentum fraction of the quark-antiquark initial state. A naive estimate allows 
us to determine that a 13 TeV pp machine with 3000 fb^ 1 of data cannot extract a statistically 
significant signal for masses, Mn r , in the range 1 — 10 TeV. However, a 100 TeV machine 
with 3000 fb _1 of integrated luminosity can yield approximately 30 events with e = 0.1 before 
a full detector simulation is implemented. In Appendix A we show an explicit analysis of the 
various contributions to the T-channel amplitude that is used to determine a series of efficient 
cuts to extract a signal for Majorana masses above 1 TeV. In particular, we show that the 
amplitudes are proportional to the jet-pair invariant mass with the strongest enhancement 
arising from contributions proportional to the Majorana mass. In addition, for momentum 
transfer, \q\ 2 , below Mjy R , the amplitudes are also proportional to the lepton pair invariant 
mass. 



M n [TeV] M n [TeV] 

(a) (b) 

Figure 10. Survival rate for a same-sign lepton final state with at least two jets for 13 and 100 TeV 
center of mass energies. 

In Figure 10 we show the performance of the detector simulation after implementing the 
default cuts introduced in Equation (3.3). For the range of masses that we are interested 
in we show the same-sign lepton reconstruction efficiency as a function of Mrt r in both 13 
and 100 TeV machines. In both cases, the overall drop in efficiency is mainly due to the 
large pseudorapidities associated with jets in the hard process since the track reconstruction 
efficiency drops significantly for values of |r/| above 2.5. In fact, from Figure 11 we can see 
how boosted jets are in comparison to the SM background. 

The T-channel dominated signal is especially interesting since in the M^ r > \q \ 2 limit, 
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Rapidity of the leading jet 


Rapidity of the subleading jet 


(a) 


(b) 


Figure 11. Pseudorapidity distribution of leading (left) and subleading (right) jets for the signal and 
SM backgrounds using a heavy Majorana neutrino with mass Mjy R = 5 TeV. 


the distribution of kinematical variables does not depend on the mass of the heavy Majorana 
neutrino. In addition, in this limit, the leptons are highly isotropic with px values signifi¬ 
cantly higher than those of the SM background. This is shown in Figure 12. Finally, the 
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0 500 1000 1500 2000 


PT of the leptons subleading lepton [GeV] 


(a) 


(b) 


Figure 12. px distribution of leading (left) and subheading (right) leptons for the signal and SM 
backgrounds using a heavy Majorana neutrino with mass M/v K = 5 TeV. 


lepton-pair separation, A Rp, and invariant mass distribution, m«, can also be used to reduce 
the background. The former is due to the fact that both leptons are back-to-back in the 
transverse plane. This is shown in Figure 13. Together with the default cuts introduced in 
Equation (3.3), we implement a series of additional cuts to extract a signal arising mainly from 
a T-channel dominated amplitude labeled T(s)13 and T(s)100 for 13 and 100 TeV machines 
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Figure 13. A Ru (left) and mu (right) distributions for the signal and SM backgrounds using a heavy 
Majorana neutrino with mass Mn r = 5 TeV. 


respectively: 


T13 

mu > 500 GeV, px.i > 200 GeV, r^| > 2, 


T100 

mu > 500 GeV, px.i > 350 GeV, ru > 2, 


Tsl3 

miihjih > 1500 GeV > e j-d > 10_5 ’ 


TslOO 

m hl 2 ,jij 2 > 2500 GeV, > 10 -5 . 

(3.11) 


While we emphasize the this strategy can yield large significances, it is very much dependent 
on the fake rate associated with the QCD background. That is, if the fake rate is higher, 
a tighter cut on the lepton px must be used. Furthermore, the lepton separation, A Ru, is 
highly correlated with mu and the lepton’s transverse momentum; we observe that it only 
reduces the signal after all other cuts are applied. 




e + e + i , i i jj 

e + p + vu jj 

p+p+j/n jj 

4 J • £ Ui 

Parton level 

cr (fb) 

0.3927 

0.7849 

0.3927 

1.79 


N 

1178 

2354 

1178 

5370 

Detector level + (3.3) 

a (fb) 

0.1187 

0.2674 

0.1187 

0.471 


N 

356 

802 

356 

1410 

T13 

a (fb) 

3.9 x 10~ B 

7.8 x 10" e 

< 10 -e 

3.0 x 10" 5 


N 

0.012 

0.003 

< 0.003 

0.089 

Tsl3 

a (fb) 

< io- b 

< 10" b 

< io- b 

2.69 x 10“ B 


N 

< 0.003 

< 0.003 

< 0.003 

0.0081 


Table 3. SM backgrounds at 13 TeV. The event number is shown at a luminosity of 3000 fb -1 . On 
the last two words we show the background reduction after applying the T13/Tsl3 cut selection in 
Equation (3.11). 


In Tables 3 and 4 we show the SM backgrounds after implementing the T(s)13 and 
T(s)100 cut selections. The Tsl3 and TslOO correspond to an additional cut on the 4-particle 
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e + e + nn jj 

e + p + vv jj 

p+p + vv jj 

4 3 ■ t'Ui 

Parton level 

cr (fb) 

3.686 

7.370 

3.694 

44.2 


N 

11058 

22110 

11082 

133000 

Detector level + (3.3) 

cr (fb) 

0.962 

2.148 

1.189 

11.6 


N 

2887 

6444 

3568 

35000 

T100 

cr (fb) 

< 10” b 

8.59 x 10” 4 

3.7 x 10” B 

7.4 x 10” 4 


N 

< 0.1 

< 0.44 

0.11 

2.58 

TslOO 

cr (fb) 

< 10” B 

< 2.0 x 10” a 

< 10” a 

2.21 x 10” 4 


N 

< 0.003 

< 0.003 

0.003 

0.663 


Table 4. SM backgrounds at 100 TeV. The event number is shown at a luminosity of 3000 fb” 1 . On 
the last two words we show the background reduction after applying the TIOO/TslOO cut selection in 
Equation (3.11). 




(a) (b) 


Figure 14. The 95% expected limit at 13 TeV and 3000 fb” 1 of integrated luminosity with varying 
fake rates on the QCD background in the same-sign electron channel (left) and same-sign muon channel 
(right). The red solid line indicates the bound from unitarity discussed in Section 2.3. 


final state invariant mass that reduces the QCD background by up to an order of magnitude 
in a 13 TeV machine. Using the cut flow discussed above, appropriate to extract a T- 
channel dominated signal; we observe that the SM irreducible background can be completely 
eliminated. In addition, implementing a fake rate of order O (lO”' 5 ) further suppresses the 
QCD background. However, the signal is proportional to the mixing angle raised to the 
fourth power, and the event yield is very small in the regions of small mixing angle. The yield 
can be enhanced with an additional cut on the 4-particle final state invariant mass by one 
order of magnitude. Using the analysis introduced in Section 3.2 we find the 95% confidence 
level fit using the x variable defined in Equation (3.9) with varying mixing angle squared, 
|U e4 | 2 , | V^ 4 1 2 . We scan over the pr of the leptons, the pseudorapidity of the subleading jet, the 
same-sign lepton invariant mass and the 4-particle invariant mass to extract the maximum 
value of x leading to a 95% expected limit on the mixing angle as a function of the Majorana 
neutrino masses. The scan is listed below: 

• !-5 < \Vj 2 \ < 3 
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(b) 


Figure 15. The 95% expected limit at 100 TeV and 3000 fb _1 of integrated luminosity with varying 
fake rates on the QCD background in the same-sign electron channel (left) and same-sign muon channel 
(right). The red solid line indicates the bound from unitarity discussed in Section 2.3. 


• 0.1 TeV < pt < 0.7 TeV 

• 0.5 TeV < mip 2 < 1. TeV 

• 0.7 TeV < mi 1 i 2 j 1 j 2 < 3.5 TeV 

The results are shown in Figure 14 for a 13 TeV collider and Figure 15 for 100 TeV. In order 
to start probing couplings allowed by unitarity, in the large Mn r region, one will have to 
wait for a 100 TeV collider and this highly depends on the ability to reduce the fake rate for 
both ej-). e and below ~ 10 -5 . For very small fake rates, a 100 TeV machine can probe 
Majorana masses up to 10 TeV with couplings as low as e e ~ 0.1 in the same-sign electron 
channel. 

3.4 Left-Right Symmetric Model: Collider reach at 13 and 100 TeV 

The presence of a charged gauged boson in the spectrum, associated with the SU(2 )r gauge 
symmetry, leads to several phenomenological differences that makes a same-sign lepton anal¬ 
ysis different from the basic Nr extension of the SM. This is due to the fact that for gauge 
boson masses, Mw R , above 1 TeV, the T-channel contribution to the production cross section 
is strongly suppressed. This makes a T-channel dominated signal difficult to extract even 
with a 100 TeV machine and 3000 fb -1 of integrated luminosity. Since we focus on the limit 
where Mr/Mn r <C 1, we can expect that BR (N —>• Wl ) —>• 0, and are not able to use the SM 
charged gauge boson mass as a means to suppress the SM backgrounds. Thus, depending on 
the mass hierarchy between the Majorana neutrino and Wr we can expect different processes 


- 22 - 












to dominate the signal: 


M Nr < m Wr 
Mn r > M Wr 


(7(pp ->• Z ± Z ± jj) dominant = & 
u{pp '? I l Jj)dominant — S 


(pp ->■ W±)iTR(W± ->• AT* 1 ± )BR(N r ->• Z± jj) 

(a) a(pp -> W±)5i2(W± ->■ Z±Z± jj) 

(b) a{pp -> W±* -> TV l ± )BR(N -»• Z± jj) 


(3.12) 


£ 


o4 

CQ 



Figure 16. Branching ratio of jj (left) and VF 1 * 1 — > jj (right) with Mw R = 2 TeV. 


The search strategy will strongly depend on whether VFr can be produced on-shell. In 
particular, a same-sign lepton plus two jets final state will be suppressed by a 4-body phase 
space factor in the region where M\y R < Mn r . We can see this in Figure 16(a) for M\y It = 2 
TeV. However, from Figure 16(b), we can see that the branching ratio of W R into two jets 
is mostly unaffected. Therefore, for Mjy R < M\y R , the pr and rj of the final state particles 
will highly depend on the pr and q of the Wr on its rest frame. However, for Mj\r R > M\y R , 
resonances can be reconstructed depending on the phase space: 

• ((I ) M\y R M t l l 2 Jij 2 ; 

• (b) M Wr = M jlj2 , M Nr = M hjjlj2 . 

To see which resonances can be used to better enhance a same-sign lepton signal, we analyze 
how (a) and (b) in Equation (3.12) behave as a function of Mn r and where the on-shell 
production of the right-handed gauge boson dominates. This is shown in Figure 17 where 
on the left column we plot the production cross section at 13 TeV as a function of Mn r for 
M\y n values of 1,2 and 5 TeV respectively while results for 100 TeV are shown on the right 
column. In the former we observe a small region for Mn r > M\y R where we cannot use the 
Ij 1 J 2 invariant mass to reconstruct the Majorana mass while in the latter, the (b) production 
mode always dominates. Thus, for a 100 TeV machine the invariant mass of the Ij i system 
can always be used to reconstruct the Majorana neutrino mass for both Mn r < M\y R and 
Mn r > M Wr . 
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Figure 17. pp —> l + l + jj cross section for varying Mn r at 13 TeV (left) and 100 TeV (right) center of 
mass energies. The top, middle and lower panels correspond to masses Mw R = 1, 2, 5 TeV respectively. 


After a fast detector simulation, most of the leptons are misidentified as jets. With 
Delphes, at 13 TeV, we find a signal survival rate of about 10 — 25%, while for a 100 TeV 
machine the survival rate is below 10% for Majorana masses above ~ 3 TeV. We note that 
this rate is unreasonably small and a more complete analysis on the simulation of different 
detector types should be used. In our study we conclude that is difficult to probe our simplified 
version of the left-right symmetric model with Mw R > 3, Mn r > 1 TeV at a 13 TeV machine 
and M\v r , Mn r > 8 TeV at 100 TeV unless we significantly reduce the SM background. In 
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what follows we only study a same-sign electron signal. A proper study of the muon chamber 
ID performance can make a same-sign muon signal just as significant and the analysis would 
proceed in the same way. 
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Figure 18. The E distribution of all final state particles with M/v B = 1 TeV (left) and M^ R = 3 
TeV. The top panel corresponds to center of mass energies of 13 TeV while the bottom panel to 100 
TeV. 


In addition to the 3- and 4-particle invariant masses discussed above, we implement an 
additional 4-particle Lorentz invariant kinematical variable that we find extremely useful to 
suppress the SM background, details which can be found in Appendix B. We label this variable 
by E and define it as 

E = evMpPfPZPgPi, (3.13) 

where -Pi, 2 , 3,4 denote the four-momentum of the different particles. The expression can be 
simplified if we define the 2 -axis to be the direction of particle p\. The effectiveness of this 
variable to reduce the SM background is depicted in Figure 18 and a qualitatively analysis is 
shown in Table 5. 

In order to enhance a left-right symmetric signal over the SM background we implement a 
series of cuts on the px of the leading lepton, the 4-particle final state invariant mass, mi 1 i 2 j 1 j 2 
and E. The cuts are chosen to maximize the statistical estimator, S/\/S + B. Furthermore, 
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E [GeV 4 ] 

e + e + vv jj (13 TeV) 

Fake (13 TeV) 

e + e + uv jj Fake (100 TeV) 

Fake (100 TeV) 

10 y 

3.367% 

0.769% 

10.954% 

3.327% 

5 x 10 9 

0.983% 

0.100% 

2.358% 

0.746% 

lCF 5 

0.350% 

0.031% 

1.058% 

0.350% 

5 x 10 10 

0.023% 

0.0027% 

0.119% 

0.062% 

lo 11- 

< 0.003% 

< 0.0002% 

0.027% 

0.023% 


Table 5. The cut efficiency for the E variable on the leading SM background for 13 and 100 TeV 
center of mass energies. 


we use mass relations between Mjy R and M\y R and two different fake rates (10 4 ,10 5 ). In 
the region where Mjy R > M\y R . the width of Nr to Wr and a lepton is given by 


r(iv -»■ iw R ) 


\Vin\ 2 M Nr 


M k 

M W R 


+ 1 


(3.14) 


and for Mjy R < M\y R they are given by 


r (Wr Jill) 
T(Wr Nl) 


A lv ,2 M w,~ M k 

te |Vwl 2 M Kr 


3-2 


M k 


M w, 


(3.15) 


Using figure 17 we see that at 13 TeV, we may focus on the region where M\y R > Mn r and 
on-shell production of the SU(2 )r gauge boson, while with a 100 TeV collider we may probe 
both mass hierarchies with on-shell Wr or Nr production. 

We calculate the collider reach in the non-degenerate region for M\y R > Mn r with 
M Nr = i M\y R (top panel) and for M\y R < Mn r with M\y R = | Mjy R (bottom panel). We 
emphasize that our choice of masses is not theoretically motivated and we use it to illustrate 
the significance of a same-sign lepton signal. The total number of events before cuts depends 
on BR(Wr -a Nl) x BR(Nr -a Ijj ) for M\y R > My R and BR(Nr -a Ijj ) for Mn r > M\y R ; 
the later which is close to 100%. These allows us to scale the cross sections accordingly except 
for when the two masses are very degenerate, since in this case the leptons in the final state 
will be very compressed. We present only the same-sign electron channel and include only the 
statistical error using the least conservative fake rate of (j^i = 10 -5 . The same-sign muon 
channel is suppressed with respect to the electron channel since within Delphes, the muon 
reconstruction efficiency is lower. A realistic muon-chamber simulation may lead to better 
results for the same-sign muon final state. Results were obtained from a scan on the py of 
the leptons, the 4-particle final state invariant mass and the E variable used to obtain the 
highest significance possible with 3000 fb _1 at center of mass energies of 13 and 100 TeV. 
Furthermore, values for V e y R , V r n r = 1 are used. The scan is listed below: 


• 10 -4 TeV 4 < E < 10 _1 TeV 4 


- 26 - 























10 7 
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(c) (d) 

Figure 19. Expected background and theory prediction with Mn r = | Mw R (top panel) at 13 TeV 
on the left and 100 TeV on the right. The lower panel corresponds to Mw R = \Mn r - 


• 0.1 TeV < px < 1- TeV 

• 0.7 TeV < mi 1 i 2 j 1 j 2 < 4.5 TeV 

It is interesting to see that even a 13 TeV collider can probe Major ana neutrino masses up to 
6 TeV and gauge boson masses up to 8 TeV. In addition, the degenerate region is also very 
interesting and we calculate the reach for My r R = M\y R — 0.1 TeV and Mn r = M\y R + 0.1 
TeV. These are shown in Figures 20 and 21 respectively. In all four cases, larger values of 
IE are required to suppress the SM background for large Wr and Nr masses. In addition, 
we observe that for a 13 TeV collider, a flat px cut at around 0.1 TeV is sufficient; while 
the 4-particle invariant mass cut increases with both the mass and the collider energy. We 
recommend that the experimental collaborations explore cuts on the px of the leading lepton 
above 0.1 TeV and on the total invariant mass between 1 and 2.5 TeV at the LHC’s run 
II, and to also consider implementing the E variable. However, we point out that we have 
assumed that the PDFs implemented in MadGraph are adequate for a 100 TeV study. An in- 
depth study is necessary, especially since the plans for the construction of a 100 TeV machine 
appears to be gaining enough momentum. 
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Figure 20. Expected background and theory prediction with Mjv r = M\y R — 0.1 TeV (top panel) at 
13 TeV on the left and 100 TeV on the right. 
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Figure 21. Expected background and theory prediction with M Nr = M Wr + 0.1 TeV (top panel) at 
13 TeV on the left and 100 TeV on the right. 


4 Discussion 

In this study we have carried out an in depth analysis of two simplified scenarios that give 
rise to a heavy right-handed Majorana neutrino. In particular we have focused on a same- 
sign lepton final state that can be probed at the LHC Run II and a future 100 TeV hadron 
collider. Our work compliments existing analyses where the major focus has been on on-shell 
production of the Majorana neutrino with masses below 1 TeV. We have extended our analysis 
to masses above 1 TeV by determining a set of optimal cuts that can be used to extract a 
T-channel dominated same-sign lepton signal. In particular, we observe that the same-sign 
lepton invariant mass is an excellent variable to suppress the SM irreducible background 
while the 4-particle final state invariant mass can be used to significantly suppress the QCD 
background using a fake rate of order O (10~ 5 ). With the set of cuts given in Equation (3.11) 
we are able to probe the same-sign electron channel beyond bounds from lepton unitarity 
with 3000 fb” 1 of data at 13 TeV. Furthermore, we observe that a significant same-sign muon 
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signal can be probe, beyond limits from lepton unitarity, with a 100 TeV machine. 

In addition, we have analyzed the sensitivity to a simplified version of the left-right 
symmetric model, where in addition to a heavy Majorana neutrino, right-handed charged 
currents mediate lepton violating interactions. Given this additional degree of freedom, we 
conduct a study based on whether Wr or Nr can be produced on-shell, and we optimize the 
collider study to extract a signal in different regions of phase space. Furthermore; our analysis 
goes beyond and compliments the experimental searches by CMS and ATLAS [46, 47] in the 
region where Mn r > 1 TeV. We find an interesting 4-particle Lorentz invariant quantity, IE, 
introduced in Equation (3.13), that can be used to suppress the background in all regions 
of phase space. This variable peaks in the presence of heavy beyond-the-SM particles. The 
analysis is carried out for a mass relation between Nr and Wr and find that in regions 
where Mjy R = Mw R /2 the LHC 13 TeV with 3000 fb -1 can reach heavy gauge boson masses 
up to ~ 8 TeV, unlike the degenerate region, where the reach is only up to 4 TeV. The 
inverse hierarchy where M^ R > M\y R is more difficult to probe since the on-shell production 
is phase space suppressed, however, the reach is up to M^ R ~ 5 TeV in the region where 
Mw r = -^Vfl/2- In addition, we have not commented on the reach of the e// jj channel, 
which unlike the single Majorana extension of the SM, is much more promising in the left- 
right symmetric model, due to the additional suppression from the ratio i n 

the ji —>• ey branching ratio given in Equation (2.18). We understand that the analysis 
in [79] shows that FCC-ee(TLEP) [80] would have a better constrain on the flavor violating 
processes through rare Z decay measurement. Being a leptonic collider it has the advantage 
of lower backgrounds but it is sensitive to lower mass heavy states. On the other hand 
hadronic colliders can probe higher masses as we have shown, but require the new couplings 
to be relatively large. Hence, this study complements the FCC-ee study of heavy neutrinos. 
Optimistically, if a signal is found for a heavy Majorana neutrino in the mass range we have 
examined it will mean that simple Type I seesaw with very small mixings is disfavored. This 
will point to either more structure in the heavy neutrino sector and/or a rich scalar sector 
which our simplified model is a truncation of. 

It is our hope that both CMS and ATLAS continue their search for heavy Majorana 
neutrinos with masses above 1 TeV early in run II. We have shown that in both simplified 
scenarios, where signal rates are small, there exist a series of interesting variables that can be 
used to probe lepton violating processes at the LHC. In addition, we encourage the collabo¬ 
rations to continue their efforts at determining the fake rate of jets into muons and electrons, 
especially at 13 TeV where the effects of pile-up will be enormous. 
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Appendices 


A T-channel amplitudes 

In this section we analyze the properties of the T-channel amplitude in the Majorana neutrino 
extension of the SM. We separate the different contributions to the cross sections by looking 
at the different contributions from partons in the initial state and jets in the final state. 
The classification in terms of flavor and color and is gauge invariant. This is shown in 
Table 6. The first two dominant contributions are o(uu —>• dd e + e + : qq — qq — d) and 


Type 

Parton Pair 

Jet Pair 

Color 

qq-qq- s 

quark-quark: same flavor 

quark-quark: same flavor 

same color 

qq-qq-d 

quark-quark: same flavor 

quark-quark: same flavor 

different color 

qQ-qq- s 

quark-quark: different flavor 

quark-quark: same flavor 

same color 

1 

1 

quark-quark: different flavor 

quark-quark: same flavor 

different color 

qq-qQ- S 

quark-quark: same flavor 

quark-quark: different flavor 

same color 

qq-qQ- S 

quark-quark: same flavor 

quark-quark: different flavor 

different color 

qQ - qQ 

quark-quark: different flavor 

quark-quark: different flavor 


qa — qa 

quark-antiquark 

quark-antiquark 


qa — qa — 0 

quark-antiquark 

quark-antiquark 

colorless 


Table 6. Classification of contributions to the T-channel mediated same-sign lepton production. 

cr(uu —>• dd e + e + : qq — qq — s), which are around the same order. The third and fourth 
dominant contributions are a(dd —» uu e + e + : qq — qq — d) and a(dd —>• uu e + e + : qq — qq — s) 
respectively. The qa — qa — 0 process can mix with the off-shell 5-channel amplitude, but 
this turns out to be subdominant. Due to the hierarchy of the CKM matrix and the richness 
of the valance quarks in the proton, the four dominant processes are sufficient to perform a 
quantitative study of the T-channel amplitude. These are shown in Figure 2. Therefore we 
can write 


a(uu —> dd e + e + : qq — qq — d) : N C (N C — 1) — | Ma ~ Mb | 2 
cr(uu -»■ dd e + e + : qq - qq - s) : N c ■ - ■ -| Ma - Mb - Me + Md| 2 , 


(A.l) 
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where N c denotes the number of colors. To simplify the analysis we make the following 
definitions: 


Ww + 2 (Piki)][ M w + 2 (P 2 fa)][M^ R + 2(p\ki) + 2(p 1 li) - 2(hki)\ ’ 

c = _zl_ 

B [Myy + 2(p 1 k 1 )][M 2 v + 2{p 2 k 2 )][M 2 NR + 2(p 1 k 1 ) + 2 ( Pl l 2 ) - 2(l 2 h)} ’ 

c , = _zi_ 

Ww + 2 (Pi k 2)][M^ + 2(p 2 k 1 )][Mff R + 2(pik 2 ) + 2(p 1 l 1 ) - 2{l l k 2 )] ’ 

c = _zl_ 

[Mfa + 2{p l k 2 )][M^ + 2(p 2 k 1 ))[M 2 NR + 2( Pl k 2 ) + 2 ( Pl l 2 ) - 2(l 2 k 2 )} ‘ 


Then the amplitude square for qq — qq — d and qq — qq — s can be written as 


\M\ 2 qq _ qq _ d = 8 M$f R g 8 (k\k 2 ) (0 Pih)( P 2l2)[C 2 A - C A C B ] - ( Pl l 2 )(h P 2 )[C 2 B - C A C B ] 

+ {P\P2){hh) ( C-A < C'B) , 

\M\ 2 qq _ qq _ s = 8 M 2 NR g\k 1 k 2 ) ((pih)(p 2 l 2 [(C A + C c ? ~ (Ca + C c )(C B + C D )] 

+ (pi/ 2 )(P2M(Cb + C D ) 2 - (Ca + C c )(C B + C D )] 

+ (P 1 P 2 )(!M(C A + Cc)(C B + Cb)D , (A.3) 

where we have used the massless approximation for jets, (k\ + k 2 ) 2 = 2{k\k 2 ). Both of the 
amplitudes are proportional to the invariant mass of the jet pair, and a stronger dependence 
comes from the Majorana mass. For large Mn r , larger than the momentum transfer, \q\ 2 , 
the equations simplify since Ca ~ C b and C c ~ Cd- In this limit the amplitudes reduce to 

\ M \Iq-qq-d = 8M N R g 8 ( k l k 2 ) (pm) {hh)^A, 

\M\ qq _ qq _ s = 8Mx R g 8 (kik 2 )(pip 2 )(lil 2 )(C A + C c ) 2 . (A.4) 

Thus, in the limit where Ca = C c, the two processes qq — qq — d and qq — qq — s yield the same 
contribution to the total cross section. In this region of parameter space, for Mjy- > \q\ 2 , 
the amplitudes are proportional to the same-sign lepton invariant mass. We make use of this 
variable to enhance a T-channel dominated signal form the SM background. 


B The E 4-particle kinematical variable 

The 4-particle final state variable introduced in Section 3.4 is defined by 

E = e^pPfPSPZPZ , (B.l) 

where Pp 2 , 3,4 denote the four-momentum of each of the particles in the final state. This 
variable can be significantly simplified in the rest frame defined by the P 2 + P 3 + Pi invariant 
mass and after assigning the z-axis to the direction of p\. In this coordinate system, the 
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4-momentum can be written as 


Pi = \Pi\ 


( 1 \ 
0 
0 

V 1 / 


Pi = \p1 I 


(l \ 

sin a cos (3 
sin a sin /3 
cos a 


P 3 = 14 * 31 


(l \ 

sin cr cos p 
sin cr sin p 
\ cos a ) 


P2+P3+P4 = 


( "* 2 , 3,4 

0 

0 


(B.2) 

where m 2 , 3,4 denotes the invariant mass of particles 2, 3 and 4, and the angles defined above 
parametrized the direction of the three-momentum with respect to the coordinate system. 
Therefore, with this choice, E can be written as 


E = -m 2 , 3 ,4 \pi | \pi | \P 3 1 sin ct sin(/3 — p) sin cr. (B.3) 

and represents the invariant mass of the particles 1 — 3 multiplied by the volume spanned by 
in three-momentum space. This variable has a large impact when the SU(2 )r gauge boson 
is on-shell. To see this, we define P\ = Pi 2 , P 2 = Pi 2 , P 3 = Pj 1 and P 4 = Pj 2 in the massless 
limit to obtain 777 ^ 2 , 3,4 = M\y R . If we choose the limit where M^ R < Mw R , then we further 
have mi 2 j 1 j 2 = M^ R . With this in mind we can write 

M jj ( m w r ~ M n r )( m n r ~ M jj) sin a sin (/^ ~ P) sin Q - 
4 [(Mjy + Mjj) + {Mjj R — M^) (cos a cos a + since sin a cos (J3 — p)) ’ 

M U M W r - M lhh^ M lhh - M h) Sin a Sin (^ ~ P) Sin 

^[( M ? 2 jij 2 + Mjj ) + ~ M jj)( cos « cos cr + sin a sin cr cos (/3 - p)) ’ 

M w R ( M hi 2 jij 2 - M N R )( M N R - M Wr) sin a sin(fl - p) sin a 
4 [( Mjj R + M^ r ) + (M^ r - M^ R )(cosacosa + sin a sin cr cos- p)) ’ 

(B-4) 

with scenarios (a) and ( b ) defined in Section 3.4. This kinematical variable provides us with a 
nice separation between the SM background and heavy beyond-the-SM particles with masses 
above a TeV. 


^ Mn r <Mw r 

IE m Nr >m Wr (a) = - 
IE m Nr >m Wr (b) = - 
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